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A theory is presented for the strong enhancement of graphene-on-substrate bandgaps by attrac- 
tive interactions mediated through phonons in a polarizable superstrate. It is demonstrated that 
gaps of up to leV can be formed for experimentally achievable values of electron-phonon coupling 
and phonon frequency. Gap enhancements range between 1 and 4, indicating possible benefits to 
graphene electronics through greater bandgap control for digital applications, lasers, LEDs and 
photovoltaics through the relatively simple application of polarizable materials such as Si02 and 

Si3N4. 

PACS numbers: 73.22.Pr 



I. INTRODUCTION 

A key goal for graphene research is the development of 
applications which require substantial bandgaps, such as 
digital transistors. Graphene monolayers have zero band 
gap, but small gaps have been observed when graphene 
is placed on substrates such as SiC (250meVjf3, gold 
on ruthenium (200meV]P and predicted for graphene on 
boron nitride (lOOmeVj^. A technique for enhancing 
these gaps up to the leV order of magnitude seen in 
silicon is crucial to the development of digital electron- 
ics. Moreover, the ability to make spatially dependent 
changes to the gap of a semiconducting material should 
open a route to new electronic devices and applications, 
and the availability of strongly tunable gaps is important 
to the development of laser diodes, LEDs, photovoltaics, 
heterojunctions and photodetectors. 

Here, I investigate gap enhancement effects due to 
interactions mediated through superstrates placed on 
graphene systems where a gap has been opened with 
a modulated potential. For example, sublattice sym- 
metry breaking has been suggested for the gap open- 
ing mechanism for the graphene on ruthenium system^. 
The specific aim is to establish if the small gap size 
of e.g. graphene on ruthenium can be increased us- 
ing attractive electron-phonon coupling to vibrations in 
a strongly polarizable superstrate. In low dimensional 
materials, strong effective electron-electron interactions 
can be induced via interaction between electrons confined 
to a plane and phonons in a neighboring layer which is 
polarizabl^. In the context of carbon systems, exper- 
iment has shown that electron-phonon interactions be- 
tween carbon nanotubes and Si02 substrates have strong 
effects on transport propertied and theory has shown 
that similar interactions account for the transport prop- 
erties of graphene on polarizable substrate^. 

Besides graphene on substrates, several alternative op- 
tions have been put forward to generate gaps in graphene. 
McCann and FalkcP^l proposed that bilayer graphene de- 
velops a gap when it is gated with an electic field, and the 
bilayer graphene gap has been observed experimentally 
by Ohta et al. using ARPEP and by Zhang et al. usmg 



infrared spcctroscopjff^. 

Electron confinement in quasi- ID structures has also 
been put forward as a solution to the generation of 
bandgaps. Graphene nanoribbons with certain edge ori- 
entations wer e theo retically hypothesised some time ago 
to posess gap j^^ l ^^ and it has been shown using ab-initio 
calculations that reduction of the nanoribbon width can 
lead to substantial gap sizes due to electron confinement, 
although for gap sizes of the order of an electronvolt, 
very narrow nanoribbons with widths of the order of 10 
A are required. Increase in nanoribbon gaps up to around 
SOOmeV due to electron confinement has been measured 
by Han et al}^ although the width variation of the rib- 
bons is large. Recent developments have allowed for the 
manufacture of high quality nanoribbons by unzipping of 
nanotubed^, and theoretical predictions have been con- 
firmed experimentally, with a measured gap of 23.8 ±3.2 
mcV for an (8,1) nanoribbon!^. 

A more extreme solution involves changing the chem- 
istry of graphene monolayers. The generation of large 
band gaps in graphene on the order of several eV has 
been hypothesize d for chemical modification with hy- 
drogen (graphane)P^E^ and fluorine (fluorographene)P3. 
While forms of graphand^ and fluorographen^^SI have 
been manufactured, and hints of bandgaps have been 
found, there is currently no consensus on the size of the 
bandgap, as it is difficult to obtain even coverage of the 
hydrogen or fluorine. Also, it is likely that the wide 
bandgaps are too large for many digital applications. Fi- 
nally, I briefly mention that graphene on some substrates 
forms Moire patterns, which also lead to a modified elec- 
tronic structur^^. 

This paper is organized as follows: a model for 
graphene sandwiched between a dielectric superstrate 
and a gap-opening substrate is introduced in Sec. [H] 
In Sec. |III| results from Hartree-Fock theory in the high 
phonon frequency limit are presented. Perturbation the- 
ory and results for low phonon frequency are presented 
in Sec. |IV[ A summary and outlook can be found in Sec. 
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FIG. 1. (Color online) Graphene-substrate-superstrate sys- 
tem annotated with interactions and sublattices. Electron- 
phonon interactions between the graphene layer and super- 
strate are poorly screened, and large interactions of strength 
fn {m) are possible. Ions in the superstrate oscillate with fre- 
quency Q. A sites have energy -|-A and B sites —A due to 
the substrate. Atoms in the graphene sheet are indexed with 
vector n and ions in the superstrate with m. 



II. MODEL 

A model Hamiltonian for the electronic properties of 
graphene on a substrate modified by the presence of a 
superstrate should have at least three components. The 
first is a mechanism for electrons to move through the 
material. The second is an electron-phonon interaction 
and the third a static potential to describe symmetry 
breaking between graphene sub-lattices by the substrate. 
The presence of superstrates also opens the possibility 
of long-range electron-phonon interactions. As such, a 
model Hamiltonian for graphene on a substrate has the 
form, 

H = -t ^ (aLcn'CT + C'^.^ana) - ^ /r^("^)?^r^£7^r,^ 
{n,n')(7 nma 

+ Y^hn{Nrr, + l/2) +Y,^r,nr,^. (1) 
m na 

The interactions between electrons in a graphene mono- 
layer and polarizable ions in the superstrate are shown 
schematically in Fig. [l] The first term in the Hamil- 
tonian describes the kinetic energy of tight binding elec- 
trons hopping in the graphene monolayer with amplitude 
t, written in momentum space as J^ki't'k'^ls'^k + <l>*kcl.af,) , 
where 0fc = — i exp(ifc.<5i) and Si are the nearest 
neighbor vectors from graphene A to B sites, = 
a{l,V3)/2, S2 = a(l,-V3)/2, ^3 = (-a,0). Electrons 
are created on graphene A sites with the operator aj^ and 
B sites with c^, and the vectors n are to atoms in the 
plane. 

The next term in the Hamiltonian describes the 
electron-phonon interaction, which has the momentum 
space form Efcq 5fcq[4-qCfc(4 + d-q) + al_^ak{bl + 

Phonons are created in the superstrate above A sites with 



and above the B sublattice with d\^, so the displace- 
ment S^rn oc {h\^-\-bm) for sitc A. The interaction strength 
fm{n) — K/[{m — n)^ + 1]'^/^ has the classic Frohhch 
form (k is a coupling constant). The effective phonon- 
mediated interaction between electrons can be charac- 
terised by the function $(n, n') = J2m fm[n]fm[n']- 
A complication of this Hamiltonian is that a basis of 
two atoms is needed to represent the honeycomb lattice, 
which is the reason for two interactions; g and g. For local 
electron-phonon coupling, g vanishes and g becomes mo- 
mentum independent. For simplicity, I study this simpli- 
fied version of the electron-phonon Hamiltonian^^, i.e the 
effective interaction is approximated as a Holstein model 
which has local interaction, <i>(n, n') = Sn.n'- The third 
term represents the energy of phonons with frequency 
r2, where Nm is the number operator for phonons. It is 
appropriate to mention the effect of the electron-phonon 
interaction on suspended graphene, which simply leads 
to renormalizations of phonon and electron modes^^ ^'^ . 

The Frohlich form for the electron-phonon interac- 
tion has been demonstrated experimentally for carbon 
nanotubes on Si02^, and it has been established to be 
the main scattering mechanism for graphene on Si02^. 
At weak electron-phonon coupling, the effects of Hol- 
stein and Frohlich interactions are qualitatively similar 
on two-dimensional lattices)^ and the Holstein form is 
used throughout this paper. It is worth noting that there 
may be quantitative changes to the results presented in 
this paper for longer range interactions. 

For completeness, I briefly discuss an alternative class 
of electron-phonon interaction where phonon motion cou- 
ples to the electron hopping, leading to a Su-Schrieffer- 
Heeger (SSH) style interaction which can lead to dimer- 
ization (and could act against the mechanism used 
here^l^Il, Xo have a strong SSH interaction, a material 
must be very flexible (e.g. the SSH model is very good for 
describing polymers). Here, since the graphene is sand- 
wiched between two other materials, the material will be 
held rigid and in-plane SSH interactions are expected to 
be much smaller than the Holstein style interaction con- 
sidered here. 

I also note that there may be some modulation of the 
strength of the electron-phonon interaction due to incom- 
mensurability of the superstrate. The effects of this in- 
commensurability are straightforward to estimate in ID, 
by computing the sum $(n,n') = Y.m fm[n]fm[n'] for 
values of n intermediate to the lattice points. It is found 
that the interaction strength varies by around ±8% of 
the average value if the superstrate is incommensurate. 
It is not expected that modulations of this magnitude 
will qualitatively change the results presented here. 

To complete the model of the graphene-substrate- 
superstrate system, the final term in the Hamiltonian 
describes interaction between electrons and a static po- 
tential. An, induced by the substrate. In particular a 
modulated potential where A sites have energy A and 
B sites —A leads to breaking of the symmetry between 
A and B sites and gives rise to a gap. In the follow- 
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FIG. 2. Schematic of the physical processes highlighted in 
the large phonon frequency limit, (a) The modulated poten- 
tial with magnitude A opens a gap. (b) The effect of Coulomb 
repulsion is to stop both electrons sitting on the same site, ef- 
fectively closing the gap. (c) The attractive phonon mediated 
electronic interaction with coupling constant A pulls electrons 
onto the same site and reduces the hopping, effectively en- 
hancing the gap. d is the difference in electron occupation 
between sites A and B. 



ing, I examine the effects of electron-phonon interaction 
on this gap. The effect of phonons on substrate induced 
gaps has not previously been studied, and as I show, the 
imphcation for the gap is significant. 



III. HIGH PHONON FREQUENCY 

To illustrate the core physical content of the model, 
I initially study the large phonon frequency limit virhere 
a Lang-Firsov canonical transformation can be used to 
derive an effective Hubbard HamiltoniarPSI 



H 



E 



'^'k'^i.Ba'^kAa 



H.c. 



(2) 

Where [/' = [/- 2tX, \e'^\^ = 3t'^ + ^'^(2 cos(fcy^/3) + 
4 cos(fcj,\/3/2) cos(3A;:r/2)) is the dispersion for the 
graphene lattice, t' — texp{—tX/V,) and A = 
$(0, 0)/2tMr2^ is the dimensionless electron-phonon cou- 
pling which is expected to be smaller than unity (M is 
the ion mass). A local Coulomb repulsion U has been 
included for completeness, although the effects of this 
in the graphene monolayer are limited to renormaliza- 
tion of the electron bands since no phase transition (to 
e.g. a Mott insulator) is measured, avoiding the need 
for more complicated treatments of the Coulomb repul- 
sion. Fig. [2] shows the basic physical processes that lead 
to gap modification in graphene on substrate systems in 
the large phonon frequency limit. Panel (a) shows the 
gap induced from the modulated potential, where elec- 
trons prefer to sit on the lower energy site B. If sufficient 
Coulomb repulsion is applied, electrons sit on different 
sites, and this leads to an effective lowering of the gap as 
shown in panel (b). Contrary to this, an electron-phonon 
interaction leads to lower energies when electrons are on 
the same site and also tends to decrease the effective hop- 
ping, which is expected to increase the gap (panel c). 

For weak coupling, the Hamiltonian may be decou- 
pled using the standard Hartree-Fock scheme, i.e. H w 
U'{niA-t)niAi + U'{niAi)niA\ - U' {niA\){niAi) ■ A mean 
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FIG. 3. Enhancement of the substrate induced graphene 
bandgap in the anti-adiabatic (large phonon frequency) limit. 
Enhancement factors of ~ 4 can be seen for medium sized 
A. Here, parameters T = 0, RSI = A = O.lt are used. 
The results here are illustrative of the core physics, and solu- 
tions of full Eliashberg style equations for physically realistic 
parameters are shown in Fig. |4] 



field solution is then taken, with {riiA) = n + d and 
{niA) = n — d. In the following, the system is half- 
filled. Minimizing the total energy with respect to d, 
a gap equation for d is obtained. 



1 



Vbz 



{Ud/2-tXd + A)/2 



A'2 



(3) 



where A' = (U - 2tX)d/2 + A is the effective bandwidth 
once interactions are taken into account, and Vbz is the 
Brillouin zone volume. This may be solved by using a 
binary search. The results for A' can be seen in Fig. 
[3] The effect of increased Coulomb repulsion is a small 
decrease in the effective gap. The effect of increased 
electron-phonon coupling A is far more dramatic, and gap 
increases are seen for all A at any U value. I note that 
the enhancement increases as the bare gap, A, decreases. 
Crucially, enhancement factors of ~ 4 can be seen for 
medium sized A which could increase the moderate gaps 
seen experimentally in graphene on substrate systems to 
the ^ leV gap sizes of silicon and germanium, with the 
caveat that changes in the form of the electron-phonon 
could quantitatively change this result. By changing the 
form of the superstrate or substrate, significant control 
could be exercised over the graphene band gap. 



IV. LOW PHONON FREQUENCY 

Given the large gap enhancement in the large phonon 
frequency limit, it is appropriate to examine if the en- 
hancement is present at low phonon frequency. It is not 
obvious that the intuitive gap enhancement seen at large 
phonon frequencies should still occur at low frequencies 
where retardation effects could lead to large pairs that 
are not localized. For low phonon frequency and weak 
coupling, low order perturbation theory can be applied. 
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I derive a set of self-consistent equations assuming the 
following ansatz for the self energy, 



iuJn{l ~ Zn) + An 

iw„(l-Z„)-A„ 



Here, the momentum independent form of the ansatz (lo- 
cal approximation) is reasonable because of localization 
by the modulated potential A and the electron-phonon 
interaction. Off diagonal terms are absent because they 
do not feature in the lowest order perturbation theory. 
The quasi-particle weight, Z„, is shorthand for Z(iuJn) 
and A„ is the gap function. Matsubara energies for 
Bosonic quantities are ujg = iTrksTs and for Fermions 
are ujn = 2-KkBT{n -\- 1/2), where T is the temperature 
and n and s are integers. 

The non-interacting graphene Green function in the 
presence of a modulated potential has the form. 



Go ^(fc,iw„) 



(4) 



The full Green function can be established using Dyson's 
equation G^"'^(fc, ia;„) = Gq ""^(/c, ia;„) — S(iw„), leading 
to. 



G ^(fc,ia;„) 



iuJnZn - A - A„ 



iio.„Zn + A + A, 



.(5) 

Substituting the expression for the Green function into 
the lowest order contribution to the self energy, 

,^^Gii{k~q,iuJn-s)do{q,uJs). 

(6) 

Here, the phonon propagator, (io(<7, Ws) = +^1)- 
The off diagonal elements of the lowest order self energy 
are zero in the case of local interaction, so they do not 
feature in the ansatz. Since there is a modulated poten- 
tial, it is necessary also to take into account the tadpole 
diagram (which has a contribution since the effects of the 
modulated potential can not simply be absorbed into the 
chemical potential) leading to the following Eliashberg 
style equations for A„ and Zn, 



An = 2t\5n - tXkeT ^ de 



Sn — ksT'^^ J de 



Dje) A'„_Miio,) 

,2 72 I A /2 I ,2 ' 

(7) 
(8) 



^2^2 A'2 -he2 

n n ' n ' ^ 



and 

_ tXksT sr-^ f D{e)uJn-sZn~sdQ{iuJs) 

^" - 1 2^ / ^^^T2 ^2 _^A'2 j_,2' 



LOn 



' uj^ Z^ + A'' + 



where Sn is the difference between the density of electrons 
on sites A and B and the full gap is A^ = A„ -I- A. Here, 



the density of states for graphene in the absence of a gap, 
D{e), has the form given in Ref. [231 Note that these gap 
equations differ from those for a superconductor. 

The equations may be solved self consistently by per- 
forming a truncated sum on Matsubara frequencies. The 
maximum Matsubara frequency was kept constant with 
the value, Wmax — 75t which is sufficiently large to 
ensure that asymptotic behavior of the gap function 
was achieved. The longest self-consistent solutions took 
around 2 weeks, increasing as 1/T and A. The following 
parameters have been chosen to match graphene with 
a modest bare gap of the order of magnitude seen in 
the systems discussed in the introduction: A = O.OSt 
corresponding to a gap of around 280meV, realistically 
achievable phonon energies of hn, = O.Olt = 28meV, 
A < 1 and temperatures on the order of room temper- 
ature ksT — O.Olt corresponding to ~ 324K. The dy- 
namical quasi-particle weight, Zn is of order unity for 
all parameters considered here and the gap equation has 
very weak frequency dependence at the values of 57 that 
have been considered. 

The gap enhancement factor A'/A is shown in Fig. 
|4] Significant gap enhancement can be seen, with a 
swift rise at intermediate A and achieving a factor 3 
at around A = 0.8. The enhancement factor increases 
slightly with decreasing A but is essentially unchanged 
by modifications to phonon frequency and temperature 
for the parameter values used here. I note that there 
may be some quantitative reduction to this enhancement 
for the longer range Frohlich interaction. Consideration 
of electron-phonon couplings of up to the order of unity 
are reasonable since even in metals such as lead, the elec- 
tron phonon coupling can be large (A ~ 1.55)2^" and in 
the absence of interplane screening, even larger couplings 
should be possible. Therefore, gap sizes relevant to digi- 
tal graphene devices should be achievable by placing ionic 
superstrates on top of e.g. the graphene on rubidium sys- 
tem, and such a system could make a good starting point 
for experimental investigations of the gap enhancement. 
It is noted that since rubidium is a conductor, an insu- 
lating material which leads to the same undressed gap 
would be necessary to make a working digital device. It 
is possible to rule out a number of materials which have 
been used for top and bottom gating in transport mea- 
surements. I note that gaps can also be induced in bi- 
layer graphen^ and graphene nanoribbon^^H^ go a similar 
mechanism of substrate mediated electron-phonon inter- 
action may enhance gaps in those systems. 



V. SUMMARY AND CONCLUSIONS 

In summary, I have presented a theory for the en- 
hancement of graphene band-gaps by polarizable super- 
strates. The theory predicts gap enhancements of up 
to 4 times from electron-electron interactions mediated 
through phonons in a polarizable ionic superstrate. The 
generation of sizable graphene bandgaps is a key problem 
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FIG. 4. (Color online) Enhancement of the substrate induced 
graphene bandgap in the adiabatic (low phonon frequency) 
limit. Gap enhancements of up to 4 are readily achievable, 
and for smaller A the gap enhancement is more pronounced 
at larger A. Realistic parameters are used: A — O.lt and 
A = O.Obt corresponding to bare band gaps of 2A = 0.56eV 
and 0.28eV respectively, t = 2.8eV, hn = O.Olt = 28meV 
and 0.02t = 56meV and ksT = 0.02t = 56meV and O.Olt = 
28meV corresponding to T = 648K and 324K and A < 1. 
In the plot fcsT = 0.02t and hQ = O.Olt unless otherwise 
specified. All energy scales are much less than the band-width 
and change in T and f2 has only a very small <1% effect on 
the gap. For such a low phonon frequency, the gap function 
A„ has a very weak frequency dependence. The frequency 
dependence gets even smaller as temperature drops. 



for the use of graphene in many technologically important 
applications. The theory shows that the relatively sim- 
ple addition of polarizable superstrates to graphene sys- 
tems can provide a way of making large enhancements to 
graphene gaps. I suggest the following recipe for experi- 
mental investigation of an enhanced graphene gap: Form 
Si02 or Si3N4 layers on top of graphene on ruthenium 
intercalated with gold. I note, however, that for digi- 
tal applications an insulating substrate is required, so if 
the effect can be verified experimentally it will be neces- 
sary to find additional substrates that cause gaps in the 
graphene spectrum. It is hoped that this work will stim- 
ulate experiment, leading to tunable gapped graphene 
systems with applications in digital electronics, LEDs, 
lasers and photovoltaics. 
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